Subwavelength imaging by a left-handed material superlens 
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In this work, a finite-difference time-domain (FDTD) method is employed to justify the super- 
lensing effect of left-handed material (LHM) slabs. Our results demonstrate that subwavelength 
resolution can be achieved by realistic LHM slabs with finite absorption and dimension. We present 
the dynamic feature of the imaging process and the dependence of physical parameters on the perfor- 
mance of the superlens. These results help to clarify the diversed FDTD results reported previously. 
We also show that the achievable resolution is limited by the absorption and thickness of the LHM 
slabs, which introduces difficulties in practical applications of the superlens. 
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I. INTRODUCTION 

In 1968, Veselago pj predicted that a planar slab of 
left-handed materials (LHMs), which possess both neg- 
ative permittivity and negative permeability, could refo- 
cus the electromagnetic (EM) waves from a point source. 
Recently, Pendry 0] extended Veslago's analysis and fur- 
ther predicted that evanescent waves, which carry sub- 
wavelength structural information of the object, could be 
amplified inside the LHM slab and reconstructed in the 
image plane without loss in amplitude. Therefore, the 
LHM slab can be used as a superlens to achieve super- 
resolution, which overcomes the diffraction limit of con- 
ventional imaging systems. However, several recent anal- 
yses 0, 0, IE IE contradicted Pendry's proposal. The 
feasibility of the superlens is still under dispute. 

There have been several attempts 0, S 13 to verify 
the superlensing effect of a LHM slab using the finite- 
difference time-domain (FDTD) method. The advantage 
of FDTD is clear. In FDTD, the results are directly ob- 
tained from Maxwell's equations and the constitutive re- 
lations of the materials, so that unnecessary assumptions 
and unessential complications can be avoided. Ziolkowski 
and Heyman Q first showed that EM waves from a line 
source could be focused paraxially by a LHM slab, but no 
stable image could be formed. They found that the image 
moved back and forth over time and sometimes vanished 
altogether. Loschialpo et al. Q observed a stable image 
of a line source formed by a LHM slab. They found that 
the image is of the order of the wavelength, showing no 
superlensing effect. However, a recent FDTD simulation 
by Cummer || reached a totally different result: sub- 
wavelength resolution of the image could be achieved by 
a LHM slab. The FDTD method is supposed to be accu- 



rate and conclusive. However, these contradictory FDTD 
results have caused more confusions. Further clarifica- 
tions are therefore needed. 

In an earlier FDTD study of ours we proposed 
a method to simulate different evanescent wave compo- 
nents separately. The method enabled us to explicitly 
study the dependence of wave behavior on different pa- 
rameters. Our simulation results provided direct numer- 
ical evidence that evanescent waves could be amplified 
in a LHM slab. The main purpose of the present work 
is to extend our previous study and investigate the con- 
tribution of evanescent waves to the image quality. We 
demonstrate that super-resolution can be achieved by a 
LHM slab with physically realizable parameters. Fur- 
thermore, we study the dependence of the image quality 
on different parameters, which may help to understand 
the limitation on the performance of the superlens. As a 
secondary motivation for our present work, we attempt 
to clarify the diversed results of previous FDTD simula- 
tions and explain why different results were obtained. 



II. NUMERICAL METHOD 

Following Ref.9, we consider a planer LHM slab which 
occupies < x < L in vacuum, with surfaces normal to 
x direction and extends to infinity in the y and z direc- 
tions. The LHM is isotropic and characterized by causal 
permittivity (e) and permeability (fi) of identical plas- 
monic form, 



e(uj) = ii(u>) = 1 



p 
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where lu p is the plasma frequency and v c is the collision 
frequency. A z-polarized electric field source is excited 
in the plane at x = —L/2 and its image is supposed to 
form in the plane of x = 3/2L. The source is sinusoidal 
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FIG. 1: Time evolution of the E z amplitude at the image 
plane for a LHM slab (L = 80 Ax w 0.14A ) with different 
values of 7, ranging from 0.005 to 0.05. The evanescent wave 
has k y = 3.0. The time evolution of the source field is also 
present for reference. 



of to w 1.1 x 10 10 rad/s and with a smooth turn-on of 30- 
periods Q • By choosing oj p and v c appropriately, we have 
e(w) = = —1 — 17 (7 denotes the loss term) in all 

our simulations. Following the method introduced in our 
previous work Q , we simulate the wave components with 
different transverse wave numbers separately by applying 
the boundary condition in the y direction for a given k v , 

E z (x,y±Ay) = E z (x,y)e* ik y A y. (2) 

Here, we restrict our simulations to evanescent waves 
with k y > fcrj (ferj = w/c), since it is the evanescent waves 
that are responsible for the subwavelength imaging. In 
the simulations, the computational space is 4000 x 1 
cells with Ax — Ay — 0.3 mm and the time step is 
At = Ax /(2c) = 0.5 ps. At the working frequency, the 
wavelength of the evanescent waves (Ao) is about 566 Ax, 
which ensures the convergence of the simulation results. 
The LHM slab is located at the center of the space in 
the x-direction. At both ends in the ^-direction, a 10- 
cell uniaxial anisotropic perfectly matched layer (UPML) 
[Tof is added to truncate the computational space. 



III. RESULTS AND DISCUSSIONS 

Upon stimulation, the source radiates an evanescent 
wave of a given k y . The evanescent wave will interact 
with the LHM slab to form an image at the image plane. 
An initial transient stage is present in the simulation, 
and only after a certain duration can a stable image with 



FIG. 2: Transfer function of a LHM slab (L = 80Ax w 
0.14Ao) for different values of 7, ranging from 0.005 to 0.05, 
as a function of normalized transverse wave number ky/ko,. 



unchanged amplitude be observed. The time required to 
reach the steady state varies dramatically, and depends 
on different simulation parameters (k y , L and 7), espe- 
cially the absorption 7. Figure 1 shows how the absorp- 
tion term 7 affects the time evolution of the image. It 
is found that the transient stage, characterized by the 
amplitude modulation, occupies longer time for smaller 
values of 7. For 7 = 0.05, it takes less than 80 peri- 
ods (T = 2ir/u>) for the EM wave to reach steady state, 
while the time required for 7 = 0.005 is about 700 pe- 
riods! Further reduction in 7 may require impractically 
long simulation time to reach the steady state. In the ex- 
treme case of zero absorption, the transient stage never 
decays and no steady state can be achieved. Our results 
are consistent with Gomez-Santos' analysis |ll|. The dy- 
namical feature of the imaging process explains why the 
FDTD simulation in Ref.3 did not show stable image. 
The authors attributed the instability of the image to 
the dispersive nature of the LHM. However, we believe 
that the unstable image comes directly from the zero or 
very small absorption used in their numerical examples. 

Next, we discuss the steady state results obtained. Fig- 
ure 2 shows the transfer function (r), defined by the ra- 
tio of the field at the image plane to that of the source, 
as a function of the normalized transverse wave number 
ky/ko. The simulation is carried out with a LHM slab of 
thickness L = 80Ax for different values of absorption. In 
the case of small absorption, for example 7 = 0.005, the 
transfer function is close to unity for small k y (|r|=0.994 
at k y — 1.5fco) and decreases slowly with increasing k y . 
Beyond a prominent peak (|r|=2.783) at k y — 4.75fco, the 
transfer function decays exponentially. The emergence of 
the peak results from the resonance of the coupled sur- 
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FIG. 3: Transfer function for LHM slabs (7 = 0.01 for all) 
of different thicknesses as a function of normalized transverse 
wave number k y /ko- 



face polaritons I9J. Due to finite absorption, the resonant 
divergence [T2I Il3| is removed. In practice, the location 
of the peak can be used as an estimate for the image res- 
olution [l4|, which is ~ Ao/5 for the case of 7 =0.005, 
showing clear evidence of the superlensing effect. The 
resonance peak is also observed for 7 =0.01 and 0.02. 
As expected, it becomes less prominent with increasing 
7. Meanwhile, the peak position shifts to smaller k y , 
though the shift is insignificant. Except for the suppres- 
sion of the resonance peak, the three curves at 7 = 0.005, 
0.01 and 0.02 show similar features and fit well with one 
another. This indicates that the superlensing effect of 
the LHM slab is insensitive to absorption change over 
certain range. However, further increase in 7 to 0.05 not 
only totally suppress the resonace peak, but also cause 
noticeable deviation in the transfer function, especially 
for large k y , showing degradation in the image resolu- 
tion. 

We also investigate the effect of thickness of the LHM 
slab (L) on the performance of the superlens. Figure 3 
shows the steady state results obtained for three LHM 
slabs of different thicknesses but with the same aborp- 
tion. The figure clearly illustrates that L has a great 



impact on the performance of the superlens. The reso- 
lution of the image can be dramatically improved when 
the thickness of the LHM slab is reduced. The image 
resolution for the LHM slabs of L = 40Ax(« 0.07A ), 
80Ax(« 0.14A ) and 160Ax(« 0.28A ) are ~ A /8, A /5 
and Ao/3, respectively. It is noted that the achieved res- 
olution is of the same order of the thickness of the slab. 
It is expected that if L is large enough (> Ao), the resolu- 
tion enhancement by the LHM slab may be totally sup- 
pressed, i.e., no evanescent waves could be reconstructed 
at the image with considerable amplitude. In fact, we 
note that this case corresponds to the numerical exam- 
ple studied in Ref.7, where L w 3.2Ao (L — 9.5cm and 
Ao = 3cm). Based on our simulation results, it is not 
surprising that the image obtained in Ref.7 showed no 
subwavelength resolution. 

The above simulation demonstrated that subwave- 
length resolution can be achieved by a realistic LHM slab, 
which is consistent with Cummer |8j. It is worth noting 
that our results are obtained directly by simulating each 
k y component, which is different from the approach used 
in Ref.8 where the amplitude of different k y components 
at the image are obtained using spatial Fourier transfor- 
mation. Our results also revealed that the resolution is 
limited by the absorption as well as the thickness of the 
slab. The absorption is unlikely to be an obstruction for 
practical application of the superlens, since we believe 
that absorption of the order of 0.01 should not be too 
difficult to achieve in artifical LHMs ^5|. However, the 
great impact of the thickness of the slab on the achiev- 
able image resolution demonstrate a stringent constraint 
on the most valuable application of the superlens in op- 
tics. For example, the superlens seems unlikely to be 
applicable in deep-submicron photolithography or nano- 
lithography 0, due to the difncluty in fabricating ex- 
tremely thin LHM slab with enough mechanical strength 
to keep the slab flat without warping. 



IV. CONCLUSION 

To conclude, we have carried out FDTD simulations to 
examine the superlensing effect of the LHM slabs. Our 
results have justified that the resolution of the image 
achieved by using the LHM slab can overcome the diffrac- 
tion limit of conventional imaging systems. However, the 
resolution of the superlens is limited by the absorption 
and the physical dimension of the LHM slabs. 
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